
SOLUTION OF DELAUNAY’S 
EQUATIONS OF TYPE H 


by 

Theodore L. Felsentreger 


April 1965 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 



SOLUTION OF DELAUNAY'S 
EQUATIONS OF TYPE H 


by 

Theodore L. Felsentreger 


SUMMARY 



An analytic solution of Delaunay's Equations of Type II 
is obtained. The results can be applied to artificial satel- 
lite problems, including satellites of small eccentricity. 
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SOLUTION OF DELAUNAY’S 
EQUATIONS OF TYPE H 

by 

Theodore L. Felsentreger 


INTRODUCTION 


Delaunay, in the expansion of his lunar theory, found it necessary to solve a particular pair of 
differential equations involving eccentricity and the argument of a trigonometric term in the dis- 
turbing function (Reference 1). These he called his "Equations of Type II"— they differ from the 
Equations of Type I in that eccentricity appears also as a divisor. They have the general form 


de 

dt 


M 1 + 


2 > 
i i 



sin 6 


d e 

dt 




(i) 


where e is the eccentricity, 6 is the argument, M is the small parameter of the problem, and N, M i? 
N i? P i are all of zero order with respect toM. Furthermore, the coefficients are assumed to be 
independent of e, 9, and time. 

These equations were also discussed by Tisserand (Reference 2) and Brown (Reference 3). 
The results derived in this paper are an extension of those presented by previous authors— the 
solutions for e cos 9 and e sin Q are complete through order e Q 4 q 3 , and the formulas for e and 9 
are also extended. In addition, the method of solution is outlined in greater detail for the conven- 
ience of the reader. 


SOLUTION FOR e cos 9, e sin e 

x 

To avoid difficulties arising from the presence of e as a divisor, we consider the variables 

x - e cos 9 , y - e sin 9 . 
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Then, from Equations (1), 


dx 

dt 


NQ 


( °° 00 \ / CO 

^~^ M i e 21 " 2 - P A e 2t “ 2 j xy - N ^1 + N i e 2i 


dy 

dt 


NQ + NQ 



e 2 


+ N 



x , 


where Q = M/N. 


A solution will be derived for x and y which is of the form 


( 2 ) 


= E n 




cos 6 


*Z> 

n = 2 


cos n 6 


e Q sin 8 


2^E„' s in n 8 [8 - 0 Q (t +c>] 


(3) 


where e 0 , c are arbitrary constants, & Q is of order zero, and E n , E n ' are of the order e 0 n q for even 
n and order e Q n q 2 for odd n. Assuming the convergence of Equations (3) is such that term-by-term 
differentiation is admissible, we have 


3F " ~ ( e o +E i) g o sing ~ /^nE n e o sinng 

n= 2 


dy 

dt 


e Q 8 o cos0 + J" 1 n E n ' 
n = 2 


0 Q cos n 6 . 


(4) 


A solution complete through order e 0 4 Q 3 necessitates the inclusion of terms with arguments 
~6, 20, 36, and 4 8 in Equations (3). Equations (3) and (4) then become 

x = E 0 + (e 0 + Ej ) cos F + E 2 cos 20 + E 3 cos 30 + E 4 cos 40 

y - e 0 sin# + Ej* sin 20 + E 3 ' sin 38 + E 4 ' sin 4 6 ( 
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and 


dx 


UX , \ — 

= - (e 0 +EJ 0 O sin 0 - 2E 2 5 0 sin 26 1 - 3E 3 e 0 sin 35 - 4E 4 0 O sin4<9 


dy _____ 

g^r = e 0 e 0 cos Q + 2E 2 6 0 cos 2 Q + 3E 3 ' & 0 cos 3# + 4E 4 ' cos 4 6 


Equations (2) are 


3^ = Nq[(Mj P ! ) xy + (M 2 -P 2 )xye 2 + (M 3 ~P 3 )xye 4 ] - N(y +Nj ye 2 +N 2 ye 4 + N 3 ye 6 ) 


g^r = NQ + NQ(Mj y 2 + Pj x 2 + M 2 y 2 e 2 + P 2 x 2 e 2 + M 3 y 2 e 4 + P 3 x 2 e 4 ) + n(x + N 2 xe 2 + N 2 xe 4 + N 3 xe 6 ) , (7) 


( 6 ) 


where 


r2 = 


(v + I e o 2 + e o E . + I E 2 2 ) + ( 2e 0 E 0 + 2E 0 El + e 0 E 2 + E, E 2 ) cos 5 
+ (T e o 2 + e o E i +2E o E 2 +e o E 3 ) cos25 + (e 0 E 2 + E 1 E 2 + 2E 0 E 3 )cos 3 5 
+ ( 2E o E 4 + e o E 3 + I E 2 2 ) cos 45 

(l e o 2 + 4 E 2 2 ) + e 0 E 2 cosd + ("T e 0 2+ e 0 E 3‘) cos25 
- e 0 E 2 cos 3 # - ^e Q E 3 + E 2 2 ) cos 4 # 


? 2 - x 2 + y 2 


xy - 


( e 0 2+E 0 2+e 0 E l + I E 2 2+ i E 2 2 ) + ( 2e 0 E O +2E 0 E l +e 0 E 2 +e 0 E 2 + Ej E a ) COS 0 
+ (e 0 E 1 + 2E 0 E 2 + e 0 E 3 + e 0 E 3 ')cos25 + (e 0 E 2 -e 0 E 2 ' + E t E 2 + 2E 0 E 3 ) cos 35 

+ ( 2 E 0 E 4 + e„ E 3 - e 0 E 3 * + \ E 2 2 - \ E 2 ' 2 ) cos 45 

( e o E o + 1 e o E 2 - T e o E 2 + T E i E 2') sin ^ 

+ (i e 0 + 1 e 0 E 1 + E 0 E 2' + \ e 0 E 3' - \ e 0 E 3 ) Sin 25 

+ [i e 0 E 2 + T e 0 E 2' + i E 1 E 2* + E 0 E 3' ) Sin ™ + ( E 0 E 4 + T e 0 E 3 + T E 2 E 2 + 1 e 0 E 3 ) Si " * 
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e “ = e o 2 ( e o 2 + 4E 0 2 + 2e o E l + E 2 + 2E o E 2 ) + 2e 0 E 0 (2e Q 2 + 2E Q 2 + 5e Q E, + 3E 0 E 2 + E 0 E 2 ' ) cos 0 

+ 2e 0 2 (E 0 2 + e Q Ej + 4E 0 E 2 ) cos 20 + 2e 0 E 0 (e Q E t + 3E 0 E 2 - E„ E 2 ' ) cos 30 

+ (2e 0 2 E 0 E 2 -2e 0 2 E 0 E 2 ')cos45 

xye 2 = | e 0 2 E 0 sin0 o (t + c) + | e 0 2 (e 0 2 +3E 0 2 + 2e 0 E 1 +6E 0 E 2 ') sin20 
+ \ e 0 3 E 0 sin 30 + (j e 0 2 E 0 E 2 + j e Q 3 E { ) sin 40 

ye 2 = e 0 (e 0 2 +E 0 2 +4 e o E i +E o E 2'" E o E 2) sin ^ 

+ (e 0 2 E 0 + e„ E„ E t + 2e Q 2 E,' + E 0 2 E‘ + e 0 E, E' + e„ E 0 E 3 ' - e Q E 0 E 3 ) sin 20 

+ 1 e o ( e o E i + 2E o E 2 + 2E 0 E 2 ')sin3^ 

+ (l e 0 2£ 2- 4 e 0 2E 2 + 4 e 0 E 1 E 2 +e 0 E 0 E 3 + T e 0 E l E 2' + E 0 E 2 E 2 + e 0 E 0 E 3 + ^ E 4 ‘ ) S in 40 
xe J = \ ( 4e 0 2 E 0 + 2E Q 3 + 6e 0 E 0 E, + e Q 2 E 2 + e Q 2 E 2 ' + 3e 0 E, E 2 + e fl E, E 2 ' + 3E 0 E 2 2 + E 0 E 2 ' 2 ) 

+ T e 0 ( 2e 0 2 + 6E 0 2 + 5e 0 E 1 + 6E 0 E 2 + 2E 0 E 2' ) COS ^ 

+ (e 0 2 E 0 + 3e 0 E 0 E, + 2e Q 2 E 2 + 3E 0 2 E, + 3e 0 E, E 2 + 3e Q E 0 E 3 + e fl E 0 E 3 ' ) cos 20 
+ J e o ( e o E i + 6 E 0 E 2 - 2E 0 E 2 ' )cos 30 + ( 3E 0 2 E 4 + 3e Q E 0 E 3 - e fl E 0 E 3 ' + \ e 0 2 E 2 
- ~2 e 0 2 E 2 ( + ~2 E o E 2 2 “ 2" E o E 2 <2 + "2 e 0 E l E 2 _ 2 e o E 1 E 2 ) cos 

X 2 e 2 =4 e o 2 ( 2e o 2 + 14E 0 2 + 7e o E 1 + 10E 0 E 2 + 4E 0 E 2 ) + \ e o E o cos F 
+ \ e o 2 ( e o 2 + 5E 0 2 + 4e 0 E 1 + 14E 0 E 2 ) cos 25 + T e o 3 E o cos # 


+ (i e 0 3 E 1 + 1 e 0 2 E 0 E 2 ” e 0 E 0 E 2 ) COS 


40 
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y2 e2 = ? e o 2 ( 2e o 2 + 2E o 2+e o E i + 4E 0 E 2 '-2E 0 E 2 ) + |e 0 3 Eo cose 

+ W (-e 0 2 - E o 2 + 2E o E 2) cos 2^ - \ e 0 3 E 0 cos 30 

~ (i e 0 3 E 1 + 1 e 0 E 0 E 2 + e 0 2 E 0 E 2 ) COS $ 

xe4 = e o 2E o( 3e o 2+6E o 2 + 9e o E i + 6 E o E 2 + 3E o E 2') +9e o 3E o 2cos5 
+ 2e o 2 E o ( e o 2 + 2E0 2 + 5e o E l + 9E o E 2 ) cos $ + e o 3 E o 2 cos $ 

+ ( e 0 3 E 0 Ej + 6e 0 2 E 0 2 E 2 - 3e 0 2 E 0 2 E 2 ' ) cos 40 

ye 4 = 3e g 3 E„ 2 sin 0 + 2e Q 2 E„ (e 0 2 + E 0 2 + 2e Q E, + 3E 0 E 2 ' ) sin 20 + e 0 3 E 0 2 sin 30 
+ ( e 0 3 E„ Ej + 3e 0 2 E 0 2 E 2 ) sin 40 

e 6 = 9e Q 4 E 0 2 + 18e Q 3 E Q 3 cos 8 + 6e Q 4 E 0 2 cos 2 0 + 2e Q 3 E 0 3 cos 36 
xye 4 = 4e Q 3 E Q 3 sin 6 + 4e Q 4 E Q 2 sin 2 8 + 2e Q 3 E Q 3 sin 30 + -j e 4 E 0 2 sin 4 8 
y 2 e 4 = y e Q 4 E Q 2 + e Q 3 E 0 3 cos 8 ~ e 4 E 0 2 cos 2 8 ~ e Q 3 E Q 3 cos 3 8 -|e 0 4 E Q 2 cos 4 8 
x 2 e 4 = e 0 4 E 0 2 + 17e 0 3 E 0 3 cos^ + 7e Q 4 E 0 2 cos 2# + 3e Q 3 E 0 3 cos 30 + y e 0 4 E 0 2 cos 40 

ye 6 = 8e Q 4 E q 3 sin 20 + e Q 4 E 0 3 sin4£ 

xe 6 = 18e Q 4 E 0 3 + 16e Q 4 E 0 3 cos 2 0 + e Q 4 E 0 3 cos 46 . 

The procedure shall be to equate the coefficients of similar trigonometric terms in Equa- 
tions (6) and (7), and then to equate similar orders of q j . 

First of all, we note from Equation (6) that the constant term of dy/dt is zero. Therefore, 
equating the constant term of dy/dt in Equation (7) to zero and grouping terms of similar order, 
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we have 


0 = (Q + E 0 ) + e 0 2(i Q M 1 + | Q P 1 + 2E 0 N 1 ) + e 0 * (-§ E 2 N, + j E,' N, + \ e 0 » QM 2 + \ e* QP 2 + 3e 0 2 E 0 N 2 ) 

+ (QE 0 2 Pj +E q 3 N 1 ) + ( e 0 QEj P, + 3e„E 0 E 1 N 1 + -j e 0 2 QE 0 2 M 2 + y e Q 2 QE 0 2 P 2 + 6 e 0 2 E 0 3 N 2 ) 

+ (-5 QEj' 2 M, + \ QE 2 2 P t + § e 0 E 1 E 2 N 1 + -j e 0 E i E 2 N i + T E o E 2 2 N i + 4 E o E a' 2 N i + T e o 3( ^i M 2 

+ e 0 2 QE 0 E 2 ‘ M 2 " \ e Q 2 QE 0 E 2 M 2 + -5 e 0 3 QEj P 2 + 4 e 0 2 Q E o E 2 P 2 + e 0 2 &o E a P 2 + 9e o E o E i N 2 


+ 6e 0 2 E 0 2 E 2 N 2 + 3e 0 2 E 0 2 E a ' N 2 + 7 e 0 4 QE 0 2 M 3 + -7 e 0 4 Q E q 2 P 3 + ™ e o* *o ” 3 ) ■ 


( 8 ) 


where the groupings are in the sequence Q, e 0 2 q, e 0 4 q, q 3 , e 0 2 q 3 , e 0 4 q 3 . To order e 0 2 q, we have 


E o = q[- 1~ l(M 1 "4N 1 +P l )e 0 2 ] 


(9) 


Equating the coefficients of cos e in the expressions for dy/dt , we obtain 


e o e o = Ne 0 (l + Nl e 0 2 ) + N(3e Q N, E 0 2 + 2e Q Q Pl E 0 + E t ) + N(e 0 QM 1 E 2 ' + e 0 QP, E 2 


+ 4 e„ 2 N, E t + 3e 0 N x E 0 E 2 + e 0 N, E fl E a * + y e Q 3 QM 2 E 0 + 7 e 3 QP 2 E 0 + 9e Q 3 N 2 E 0 2 ) . 


( 10 ) 


Equating the coefficients of sin e in the formulas for dx/dt yields 


“( e o +E i)«o = ^ Ne 0 (l +Nl e 0 2 ) + N [e 0 Q(M 1 -P 1 )E 0 - e oNl E 0 2 ] 


+ N[4e 0 Q(M 1 -P 1 )E 2 '- 4 e o Q(Mi" p i) E 2‘ 4 e Q 2 n i E i ” e o n i E o E 2 +e o N i E o E 2 


+ T e o 3 Q(M2- p 2) E o-3e 0 3 N 2 E 0 2 ] • 


(ID 
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In the last two equations, the groupings are in the order e 0 , e 0 3 , e 0 q 2 , e 0 3 q 2 . To obtain approxi- 
mations to order e 0 Q 2 for E, and Q 2 for 0 O , set 

e o*o = Ne 0 (l +Nl e 0 2 ) + N(3e oNl E 0 2 + 2e 0 QP 1 £„+£,) 

-( e o +E i) 0 o = " N e o(l +N i e o 2 ) + N [ e oQ( M i- p i)E„-e 0 N lEo 2 ] , 

obtaining 

El - "5" Q 2 (Mj _ 2N 1 +P 1 )e 0 (12) 

= N[(l+N 1 e 0 2 ) + 4 Q 2 (M 1+ 4N 1 -3P 1 )] • (13) 

Now, setting the coefficients of sin 26 in the equations for dx/dt equal results in the expression 


2E 2 0 o = N[|e 0 2 Q(M 1 -P 1 )-e 0 2 N 1 E 0 -E 2 '] + n[j e 0 4 Q(M 2 - P 2 ) - 2e 0 2 N, E 2 ' - 2e 0 « N 2 E„] 

+ N [<K M l - P i ) E 0 E 2 ' + y e„ Q(Mj - P,) E, - e 0 E 0 E 2 - N, E q 2 E 2 ' + | e Q 2 Q(M 2 -P 2 ) E 0 2 - 2e 0 2 N 2 E 0 3 ] 


+ "B e o Q( M i - P i) E 3 ' - \ e 0 QK " P t ) E 3 " e o Ej E 2 ' - e„ N, E„ E 3 ' + e, 3 q(M 2 - P 2 ) Ej 


+ 3e 0 2 Q(M 2 -P 2 )E 0 E 2 ' -4e 0 3 N 2 E oEl -6e 0 2 N 2 E 0 2 E 2 + e 0 Nj E Q E 3 + 4e 0 < Q(M3 - P 3 ) E 0 2 - 8e Q < N 3 E„ 3 ] . (14) 

Equating the coefficients of cos 20 in the equations for dy/dt, we get 

“2 B o = n[- | e 0 2 Q(Mj - Pj + e 0 2 Nj E 0 + E 2 ] + n[- \ e„< Q (m 2 - P 2 ) + 2e Q 2 N, E, + 2e^ N 2 E„ ] 

+ N b&i E o E 2 + e „ CP ! E 2 + 3e 0 N, E 0 E, + 3N, E 0 2 E 2 - \ e 0 2 Q (m 2 - SP 2 ) E„ 2 + 4e 0 2 N 2 E„ 3 ] 

+ N[e fl QMj E 3 ' + e 0 QP, E, + 3e oNl Ej E 2 + 3e 0 N 3 E 0 E, + e Q N, E 0 E a ' + 2e 0 3 QP.Ej + e 0 2 Q{M 2 + 7P 2 )E 0 E 2 


+ 10e 0 3 N 2 E 0 E, + 18e 0 2 N, E 0 2 E 2 - e * QM 3 E 0 2 + 7e* QP 3 E 0 2 + 16e Q « N 3 E 0 3 ] - 


(15) 
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In Equations (14) and (15), the orders are e 0 2 Q, e 0 4 Q , e Q 2 Q 3 , and e 0 4 Q 3 . Making use of Equations (9) 
and (13), to order e 0 2 Q the Equations (14) and (15) are 


2E 2 + E 2 ' 


= \ e 0 2 Q(M 1 -P 1 ) + e 0 2 QN l 


E 2 - 2E 2 ' 


= '|e 0 2 Q(M 1 -P 1 ) + e^QN, , 


the solution being 


E 2 = E 2 ' = - ^Q(M 1 + 2N 1 -P 1 )e„ 2 . 


(16) 


For the coefficients of sin 39 and cos 36 (order e 0 3 Q 2 ), we have, respectively, 

-3E 3 0 o = NUe 0 Q(M 1 -P 1 )(E 2 +E 2 ')-^e 0 N 1 (e 0 E 1 + 2E 0 E 2 + 2E 0 E 2 ') + ^e 0 3 Q(M 2 -P 2 )E 0 -e 0 3 N 2 E 0 2 -E 3 '] (17) 


3E 3 ' = n[- e 0 QM, E 2 + e 0 QP, E 2 + - j e Q N, (e 0 E, + 6E 0 E 2 - 2E 0 E 2 ' ) - \ e 0 3 q(m 2 - P 2 ) E q + e Q 3 N 2 E„ 2 + E 3 ] . 


(18) 


Employing Equations (9), (12), (13), and (16), we deduce 


- 3E 3 + E 3 ' 


E 3 " 3E 3 


= - j e 0 3 Q 2 (2M, 2 +6Nj 2 + 2P, 2 + 9M 2 N, - 4M t P 3 - 7NjP, + 2M 2 + 4N 2 - 2P 2 ) 


which yields 


E 3 “ E 3' 


= Q 2 (2Mj 2 + 6N 2 + 2P, 2 + 9M,Nj - 4M, Pj - 7 NjPj + 2M 2 + 4N 2 - 2P 2 )e Q 3 . (19) 


Finally, we obtain the following equations for E 4 and E 4 ' (orders e 0 4 Q and e Q 4 q 3 ): 


4 E 4 0 O 


n(e 4 ‘ + \ e o N, E 2 - \ e 0 2 N, E 2 ' ) + n[q(m x " P x ) (e o E 4 + j e Q E,' + J % E 3 + \ E 2 E 2 ') 


Q(M 2 "P 2 )(f e 0 2E o£ 2+ ie 0 3 E 1 )-N 1 (|e 0 E 1 E 2 +e 0 E 0 E 3 + 4 e o E l E 2 + E o E 2 E 2 


+ e o E o E 3 + E o 2 E V ) - N 2 ( e o 3 E o E , + 3e 0 2 E 0 2 E 2 ) + \ Q(M 3 - P 3 ) e fl 4 E 0 2 - N 3 e Q 4 E Q 3 _ 


( 20 ) 
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4E; 0 o - n(e 4 + \ e 0 2 N, E 2 - | e 0 2 N, E 2 ') + n[- QM, (e # E s ' + | E 2 ' 2 ) + (J» l (2E 0 E 4 + e fl E 3 + | E 2 2 ) 

- QM 2 (I e 0 3 E 1 + T e 0 2 E 0 E 2 + e 0 2 E 0 E 2') + <?*2 (i e Q 3 E 1 + 4 e 0 2 E 0 E 2 " e 0 2 E 0 E 2 ') + N ! (^0 E 4 


+ 3e 0 E 0 E 3 -e 0 E 0 E 3 ' +|e 0 E 2 2 - 4 E o E 2 2 + I e o E i E 2 " 4 e o E i E 2 ') + N 2 ( e 0 3 E o E i 


+ 6e 0 2 E 0 2 E 2 - 3e 0 2 E 0 2 E 2 ') “ \ QM 3 e 0 4 E 0 2 + \ QP 3 e„ 4 E 0 2 + N 3 e* E 0 3 ] . 

From these equations, it is concluded that the e 0 4 Q parts of E 4 and E 4 are both zero. 

Now, let 

£E 0 = e 0 4 QA + Q 3 B + e 0 2 Q 3 C . 

Substituting the expressions in Equations (9), (12), (13), and (16) into Equation (8), we get 
O = e 0 4 q[a ~ 4 N l ( 3M 1 - 6 Nj +Pj ) + § (M 2 + P 2 ) “ 3N 2 ] + Q 3 (b + Pj - Nj) 

+ e 0 2 Q 3 [c + 2N j B + 9N X 2 + f Pj 2 -3Mj N, + -§11^-0^ + \ M 2 - 6N 2 + \ P 2 ] 
from which A, B, and C can be determined. The result is 
££ 0 = 4Q(3M 1 N 1 -6N 2 +N 1 P 1 -M 2 + 6N 2 -P 2 )e 0 4 + Q 3 (N 1 -P 1 ) 


( 21 ) 


+ ^Q 3 (-22N I 2 -a > 1 2 +6M 1 N,-3J I P 1 +20N 1 P,-M 2 + 12N 2 -7P 2 )e 0 2 . (22) 

Let A0 O be the e 0 2 Q 2 part of d Q and AE, the e 0 3 Q 2 part of E,. With the aid of Equations (9), (12), 
(13), and (16), Equations (10) and (11) become 

e o A0 o - NAEj = \ e 0 3 Q 2 N(-2M 2 -42N 1 2 -2P/+21M 1 N,-4M 1 P 1 + 21N, P 3 - 2M 2 + 36N 2 - 14P 2 ) 

-e o A0 o -NAEj = ^-e Q 3 Q 2 n(- a * 2 + 14N, 2 + 2P 2 + 5Mj N t - IINjPj -6M 2 - 1^J 2 + 6P 2 ) , 
from which it follows that 

A0 O = ^ Q 2 N(- 28N, 2 - 2Pj 2 +8Mj Nj - avijPj + 16Nj Pj+ 2M 2 + 24N 2 - 10P 2 )e 0 2 (23) 
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AEj = | Q 2 (2Mj 2 + 14N, 2 - 13MjN, + 2M! Pj - 5 NjPj +4M 2 - 12N 2 +4P 2 ) e 0 3 . 


( 24 ) 


Next, let 

AE 2 = e 0 4 QD + e 0 2 Q 3 F , AE 2 ' = e Q 4 QG + e 0 2 Q 3 H • 

Substituting the approximations obtained so far for 6 > 0 , E 0 , e 3 , and E 2 into Equations ( 14 ) and ( 15 ), 
there results 

e 0 4 Q(- 2D + G) + e 0 2 Q 3 (- 2F + H) = ^ e Q 4 Qc 4N t 2 + Mj N, + N t Pj + M 2 + 4N 2 - P 2 ) 

+ \ e 0 2 Q 3 (M, 2 - 20N 2 - 5Pj 2 - 6M t N t + 4M 2 P 3 + 22N t P t + 6M 2 + 8N 2 - 6P 2 ) 


e 0 4 Q(-D+ 2G) + e 0 2 Q 3 (-F+ 2H) = \ e 0 4 Q(4N 2 -Mj N, - N, P t -M 2 - 4N 2 +P 2 ) 

+ \ e 0 2 Q 3 (M, 2 +10N! 2 + 2P 3 2 -MjPj - 10Nj P 1 ~ M 2 - 8N 2 + 5P 2 ) , 

from which D, F, G, and H can be computed. Thus, 

£E 2 = ^Q('W 2 -M 1 N 1 -N 1 P 1 -M 2 -^l 2 +P 2 )e 4 +|Q 3 (3CN 1 2 + 7P 1 2 +6M 1 N 1 -SM 1 P 1 -3^ 1 P 1 -7M 2 -16N 2 + llP 2 )e 0 2 (25 ) 


Z£ 2 ' = ^Q(^ 1 2 -M 1 N 1 -N 1 P,-M 2 -^ 2 +P 2 )e 0 4 + + 60N 2 + 13> 2 +(M 1 N, ~ 8MjP, - 62NjPj - 1CM 2 - 4CN 2 + 26P 2 )e 0 2 . (26) 


Finally, the e 0 4 Q 3 parts of E 0 , E 2 , e 2 ' , e 4 , and e 4 ' are obtained by substituting into Equations (8), 
( 14 ), ( 15 ), ( 20 ), and ( 21 ) the results gotten so far. The complete expressions for E 0 , & 0 , etc. are as 
follows: 


E o = q[" 1"I (Mi-4N 1 +P 1 ) e 0 2 + i (sM^-GN^+NjP^ 


M 2 + 6N 2 


-P 2 ) e o 4 ] 


+ Q 3 [(Nj - Pj) + -J (- 22Nj 2 - 3Pj 2 + 6M t Nj - 3M t Pj + 20N 3 P t - M 2 + 12N 2 - 7P 2 ) e Q 2 - (M t 3 - 436Nj 3 


+ 3P X 3 - 27Mj 2 Nj + 5M t 2 P 3 + 240M t N 2 + 7Mj P t 2 + 356Nj 2 P 3 - 65Nj P, 2 - 124M, Nj Pj + 7Mj M 2 


- 144MjN 2 + 49MjP 2 - 66NjM 2 + 528N,N 2 


- 174Nj P 2 + 19P t M 2 - 168Pj N 2 + 37P t P 2 + 12M 3 - 144N 3 


+ 6° P 3) e o 4 J 


10 



<9 0 - N (l + Nj e 0 2 ) + NQ 2 [y (M, + 4N, - 3Pj) + ^ (" l^! 2 -P, 2 +4 MjN! -M,P j + 8NjP, +M 2 + 


12N 2 -5P 2 )e 0 2 ] 


E i = Q 2 [i(Mi-2N,+P 1 )e 0 + | (2M 1 2 +14N 1 2 -13M 1 N 1 +2M 1 P 1 -5N,P 1 +4M 2 -12N 2 +4P 2 )e 0 3 ] 


e 2 □ 


q[- T ( M i + - p i) e o 2 + T («i 2 " M i N i - N i P, - M 2 - 4N 2 + p 2 ) e Q 4 ] 

+ 0 3 [^ (30N, 2 + 7Pj 2 + 6M, Nj -5Mj P, - 32N, Pj - 7M 2 - 16N 2 + llP 2 )e 0 2 - jj (l8M, 3 + 2268N, 3 - 18P*. 


F ' = 

"2 


- 93Mj Nj + 6M/ Pj “ 624Mj N 2 2 - 30M 2 P 2 - 2060N/ P 2 + 433N X P 2 2 + 388M 1 N 2 P, + 138M 2 M 2 
+ 300Mj N 2 - 114Mj P 2 - 256Nj M 2 ~ 2920N X N 2 + 1088N 2 P 2 + 90P 2 M 2 + 1092P 2 N 2 - 306P 2 P 2 

+ 216M 3 + 768N 3 - 360P 3 ) e Q 4 ] 

Q[- \ ( M i + »i - ] Pi ) e o 2 + J (^i 2 '“iNr P, - M 2 - 4N 2 + P 2 ) e„ 4 ] 

+ Q 3 \j2 (S*, 2 + 60N, 2 + 13P, 2 + 6Mj N, - 8Mj Pj - 62N, Pj - 10M 2 - 40N 2 + 26P 2 ) e Q 2 - (2268Nj 3 - 18Pj 3 

+ 39M j 2 Nj - 6Mj 2 Pj - 780Mj Nj 2 - 24M } Pj 2 - 2008N j 2 Pj + 413Nj Pj 2 + 404Mj N, Pj + 24Mj M 2 


+ 564Mj N 2 - 192M, P 2 - 104Nj M 2 - 3128Nj N 2 + 1144Nj P 2 + 72Pj M 2 + 1068Pj N 2 -288Pj P 2 

+ 144M 3 + 960N 3 - 432P 3 ) e Q 4 J 


E 3 = E 3 


= ^ Q 2 (2M t 2 + 6Nj 2 + 2Pj 2 + 9 MjNj -4MjPj- 7 NjPj + 2M 2 + 4N 2 - 2P 2 ) e Q 3 


e 4 = e; 


- 72 Q 3 (9Mj 3 + 36N j 3 - 9P , 3 + 66 Mj 2 Nj - 27Mj 2 Pj + 102M t Nj 2 + 27M, Pj 2 - 70N, 2 Pj + 44Nj Pj 2 


- lOaij N,Pj + 27M, Mj + 6CM, N 2 - 27Mj P 2 + 40J, M 2 + 64Sj N 2 - 38Nj P 2 - 21?,!^ - Sffj N 2 + 21P, P 2 + 13J 3 (27) 

+ 2^ 3 -iy 3 )e 0 4 . 

SOLUTION FOR e, 6 
Now, let 


Q = 6 + Qa j sin 6 + Q 2 a 2 sin 20 + Qa 3 sin 36 - 6 + 86 


( 28 ) 
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e 


(29) 


- e ! + Q^i cos ^ + Q 2 /3 2 cos 2# + 0^3 cos 3# , 

where 



We will determine a. and /3. (i = i, 2, 3) by forming expressions for e cos e and e sin# from 
Equations (28) and (29), and then use Equations (27). 

First of all, keeping terms through orders Q 3 and l/e 0 4 , form 


sin# - sin 6 cos S# + cos 6 sin S# 

~ sin # [l - ( S#) 2 ] + cos 6 |s# - -g- ( S#) 3 ] 

- ^1 - Q 2 CLj 2 + ^ Q 2 a 2 + \ Q 2 CL 1 a 3 ) sin # + (y Q a i ~ J4 Q 3 a i 3 “ \ Q 3 a i a 2 + T Q 0 ^ ) sin 2 ^ 

+ Q 2 a i + \ Q 2 a 2 " \ Q 2 a i a 3) sin ^ 

cos 9 = cos # cos 8# - sin # sin S# 

~ cos # [l - ( S #) 2 J - sin 0 j§# - ( §#) 3 ] 

= ( _ ^ Q 3a , 3 " i Q 3a i a 2 ) + i <3 2a 2 " ?Q 2 a i a 3 ) cos e 

+ llQ“r 12 Q 3 ai 3 -T^ 3 ) cos 2F + (^QW + TQ 2a 2- 

Then, retaining terms through orders Q 3 and l/e 0 3 , 


30 


e cos 6 ~ 


[ e i (~ 2 + 16 Q 3 a l 3 " T Q 3 a i a 2 ) + I 0^1 l 1 " I Q 2 a i 2 ” 4 Q 2 a s) + {q 3 a !^ 2 ] 

+ [ e i ( 1_ -g-Q 2 a ! 2 " I Q 2a 2 ■ ? Q 2 “! a 3 ) 'T ( 5 2a i/ 3 1 + I ( 5 2 ^ + T ( 3 2a i^' TQ 2 *3/^] cos 
+ je, (j Qa J - yQa 3 - ^ Q 3 a t 3 ) + \ - 1 Q 3 “i ^2 + I 25 

+ [ e i(l Q 2 a i 2 + 5Q 2a 2 _ {q 2 a i a 3 ) + | QSid Q“i 'lQ a 3 ) +^Q 2 ^ 2 -^Q^j/SjJcos 30 (30) 
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e sin Q - (l - Q 2 ai + J Q 2 a 2 + ^ Q 2 a 1 a 3 ) + 1 Q“x + T ^3) ' £3] sin 6 

+ [ej(-| Qaj + \ Qa 3 - 33 Q 3 I 1 ? 3 a i a 2 ) + T Q 5 ! I 1 " \ QW +Q 2 cl 2 ) ~ \ Qfl 3 Jsin 20 

+ [ e i(F Q 2 a l 2 + \ Q 2 a 2 " \ Q 2a l a s) + Q 0- 3 ) ^" J Q 2 /^] sin 30 ■ (31) 

Comparing Equations (30) and (31) with Equations (5) and grouping terms of the same order, 
we obtain 


E 0 2 e x a i) + Q 3 (l6 e i a i 3 4 e i a i a 2 16 a l 4 a 2 P\ + 4 a i ^2 ) 


(32) 


; o + E i = e l + Q 2 (“¥ e l a l 2 _ T e 1 a 2 -|a 1 / 3 1 + |/ 3 2 )+ Q 2 (- ■}e 1 a 1 a 3 + J (33) 


; 1 + Q 2 (" I e i a i 2+ I e i a 2 + 4 a i^i‘ T^ 2 ) + Q 2 (i e i a i a 3 + 4 a 3 ^i 4 a i' e s) (34) 


E 2 = lQ(^i +e x a i) + Q 3 (“ U e i a i 3 ~ i e i a 3“T a x^2 + i^a) 


(35) 


E 2' 


- TQ(^I +e i a x) + Q 3 (- l4 e i a i~i e x a x a 2 + T e i a 3~ S a X 2 ^I + i a 2' e X _ 1 ^ 3 ) 


(36) 


E 3 - q 2 (I e x a x 2 + 1 e x a 2 + i a x^x + i ^ 2 ) + Q 2 ( _ T e x a i a 3 - \ a 3 ?r h a x^ 3 ) 


(37) 


E 3 ' < ? 2 (l e i a i 2 + i e x a 2 + i a i/ 3 I + '5/ 3 2) + Q 2 ( _ i e x a x a 3 + T a 3' 6 l) 


Equation (34) yields a first approximation to 


Equations (32) and (35) give first approximations to a 1 and / 3 1 : 


Pi = F)( E 2 +E o) = - 1 - (Mx-Nx)e 0 2 - (N^-MjN^M.-Nje^ 


(38) 


(39) 


a x = ^( E 2 - E o) = T; [l-(3N 1 -P 1 )e 0 2 + (5Ni 2 -2MxN I -N 1 P 1 -5N 2 +P2)e 0 ''] (40) 


(41) 
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Letting Ae t equal the q 2 part of e^from Equations (33), (34), (37), (39), (40), and (41) we have 

-I e o a 2 + I' S 2 = gii^ [-3 + 2(5N 1 "P 1 )e 0 2 + (2M 2 -23N l 2 -P 2 +3M 1 N 1 + 7N 1 P 1 +2M 2 + 10N 2 )e;] (42) 

4e 0 a 2 + -j/3 2 ^ g§^ [l + 2^ - N t ) e Q 2 + (2M t 2 + 5N t 2 + P, 2 + M, N, - 2M, P, - 3N, P, + 4M 2 + 2N 2 - 2P 2 ) e*] (43) 

Ae i = Q 2 (f e o a i 2 " 7 a i^i~T e o a 2 + ' ( 44 ) 

The solution of Equations (42) and (43) is 

*2 = A [2+(Mi"6N l +P 1 )e 0 2 t (l4N 1 2 +P 1 2 -M 1 N 1 -M 1 P l -5N I P 1 +M 2 -4N 2 -P 2 )e;] (45) 

He o 

P 2 = 4 ^ [- 1 + (M 1 + 4N,-P l ) e 0 2 + (2M 1 2 -9N 1 2 + 2M 1 N 1 -M,P 1 + 2N l P 1 + 3M 2 + 6t< 2 -P 2 )e 0 ‘’] . (46) 

Substituting these results into Equation (44) gives 

Ae i = ^[ 2 + 2 (M,- 8Nj + 3Pj) e 0 2 + (2Mj 2 + 52N X 2 + 2Pj 2 - 21M,N, + 2M, P t - 21N t Pj + 4M 2 - 32N 2 + 8P 2 )e 0 4 ].(47) 
Now, from the Q 2 / e o parts of Equations (37) and (38), we see that 

@3 ~ e 0 a 3 ' 

Therefore, letting Aa 1( A S l designate the Q 2 parts of a l and p lf Equations (32), (35), and (36) yield 

5Q(A3,-‘oto,) = [l-(2M 1 +N 1 +9P 1 )e 0 2 -(M 1 2 + 88N l 2 + 18P 1 2 -35M 1 N 1 + 22M 1 P 1 

- 114NJPJ +10M 2 -61N 2 +47P 2 )e 0 4 J (48) 

hK +e 0 Aa 1 )+ Q 3 / 3 3 - ^ [-8 + 24(3N,-P 1 ) e 0 2 + (6M 2 -78N 1 2 + 38P 2 + 123M 1 N 1 

- 52Mj Pj - 109Nj Pj - 50M 2 + 4N 2 + 58P 2 ) e* J (49) 
y Q(AAj + e fl to, )~ Q 3 fi 3 ° -^j[l4 + 6(2M 1 -5N 1 -P 1 )e 0 2 + (l2M 1 2 + 222N 1 2 + 46P 1 2 + 63M l N 1 

- 56M,P, - 191N, P x - 34M 2 - 118N 2 + 68P 2 )e 0 4 J (50) 
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The solution of these three equations is 

Aaj = ^ [(M, + 2N I +Pi) + (Mj 2 +28N* +8P, 2 -M 1 N, + M,P 1 - 41N,P,-M 2 -20N 2 + 17P 2 ) e 0 2 ] (51) 

AS. = -^•ri+(3N 1 -7P 1 )e 0 2 + (M 2 -32N 2 -2P 2 +33M 1 N 1 -20M 1 P 1 +32N 1 P 1 -12M 2 + 21N 2 -13P 2 )e 0 '*] (52) 

8e 0 L 

jS 3 = [-11-3(2M 1 -17N, + 3P,) e 0 2 -(3M 1 2 + 150N 2 +4P, 2 -30M 1 N 1 - 2M, P, - 41N, P, 

+ 8M 2 ~61N 2 + 5P 2 )e 0 4 ] (53) 

a 3 ~ ~ 'e^ j @3 ' ( 54 ) 

The complete expressions for e, and the a., /3 i are as follows: 

q 2 r- 

e i = e 0 [2 + 2(11, -8N 1 +3P,)e ( 2 + (2M 1 2 + 52N 1 2 + 2P, 2 -21M 1 N 1 + 2»» 1 P l -21N 1 P 1 

+ 4M 2 - 32N 2 + 8P 2 )e 0 4 ] 

a, - ^-[l"(3N 1 -P 1 )e 0 2 + (5N 1 2 -2M 1 N 1 -N 1 P 1 -SN 2 +P 2 )e 0 4 ] 

+ ^ [(M, + 2N, + P, ) + CM, 2 + 28N 2 + 8P, 2 - M, N, + M, P, - 41N, P, - M 2 - 20N 2 + 17P 2 ) e Q 2 ] 

Pi = - 1 ' (**i-Ni)e 0 2 - (N l 2 -M 1 N 1 +M I -N a )e« + [l + (3N, - 7P,) e 0 2 

+ (M 2 - 32N, 2 - 2P, 2 + 33M, N, - 20M, P, + 32N, P, - 12M 2 + 21N 2 - 13P 2 ) e 0 4 ] 

a 3 = — L T |’-U-3(2l»i-17N ) +3P 1 )e o 2 - v 3M 1 2 +150N 2 +4P 2 -30»l 1 Ni-ai,Pi-41NiPi+8ll 2 
4» e 0 

-61N 2 + SP 2 )e 0 4 ] 

ft 3 = - e 0 a 3 ' ( 55 ) 
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